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Abstract 

We study a strongly coupled system consisting of a parabolic equation and a singular Hamilton-Jacobi equation in one 
space dimension. This system describes the dynamics of dislocation densities in a material submitted to an exterior 
applied stress. Our system is a natural extension of that studied in [15] where the applied stress was set to be zero. 
The equations are written on a bounded interval with Dirichlet boundary conditions and require special attention to 
the boundary. We prove a result of global existence of a solution. The method of the proof consists in considering 
first a parabolic regularization of the full system, and then passing to the limit. For this regularized system, a result 
of global existence and uniqueness of a solution has been given in [16]. We show some uniform bounds on this solution 
which uses in particular an entropy estimate for the densities. 
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1 Introduction 

1.1 Physical motivation and setting of the problem 

In [12], Groma, Czikor and Zaiser have proposed a model describing the dynamics of dislocation densities. 
Dislocations are defects in crystals that move when a stress field is applied on the material. These defects are 
one of the main explanations of the elastoviscoplasticity behavior of metals (see [8] and [9] for various models 
relating dislocations and elastoviscoplastic properties of metals). This model has been introduced in order to 
describe the possible accumulation of dislocations on the boundary layer of a bounded channel. Dislocations 
are distinguished by the sign of their Burgers vector ±b (see [13] for a description of the Burgers vector). 
More precisely, let us call 9^ and 9~ , the densities of the positive and negative dislocations respectively. For 

xel := (-1,1), 
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and t € (0,T), for some time T > 0, the non-negative quantities 0^{x,t) and {x,t) are governed by the 
following system (see [12]): 



e+ + e- 



0+ + e- 



— T 



in 7x(0,T), 
in /x(0,T), 



(1.1) 



where r is the applied stress field which is assumed to be constant. Here the term Tb = l+,n- is called the 



e++e- 



back stress and can be interpreted as the contribution to the stress of the short-range interactions between 
dislocations. If r = and 9~ = (resp. 9^ = 0) then 0+ (resp. 6~) satisfies the usual heat equation. More 
generally, the back stress is proportional to the gradient of the effective dislocation density — 9~ , with a 
diffusion coefficient which is g+^g- ■ In fact, system (1.1) is a model for a 2D channel with coordinates {x,y) 



that is invariant in the y-direction (see Figure 1). The channel is bounded by walls that are impenetrable 
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Figure 1: Geometry of the crystal. 



by dislocations (i.e., the plastic deformation in the walls is zero), 
represented by the zero flux condition, i.e. 



In this case the boundary conditions are 



9+ + 9- 



— T = 0, at 



±1. 



(1.2) 



For related literature, let us mention the work of Groma-Balogh [11], where the back stress was neglected. For 
the model described in [11], we refer the reader to [5, 6] for a one-dimensional mathematical and numerical 
study, and to [4] for a two-dimensional existence result. The special case r = for system (1.1) has been 
studied in [15], where a result of existence and uniqueness has been proved. In the present paper we study 
the case for general constant r. 



1.2 Setting of the problem 

We consider an integrated form of (1.1) and we let 

Px=0''^, P = P^-P~ and K = + 
to obtain, for special values of the constants of integration, the following system in terms of p and k: 



KtKx = PtPx 
Pt — Pxx ~ Tl^x 



on /x(0,T) 
on /x(0,T), 



(1.3) 



2 



with the initial conditions: 

k{x,0)=k°{x) and p(a;, 0) = p°(a;). (1.4) 

To formulate heuristically the boundary conditions at the walls located at a; = ±1, we first suppose that 
K^; ^ at a; = ±1. We recall that the dislocation fluxes at the walls must be zero, which require (1.2). 

Rewriting system (1.3) in terms of p, k and $ = 0+_^_0- — t, we get 



pt = Kx^. 



(1.5) 



From (1.2) and (1.5), we deduce that 

Pi(.T, .) = Kt(.T, .) = for x^±l. (1.6) 
In this case, we consider the following boundary conditions: 



k(x, .) = X for X = ±1 
p{x, .) = for x = ±1, 



(1.7) 



where we have taken the zero normalization for p on the boundary of the interval. 

The non-negativity oi 6^ >0 reduces in terms of p and k to the following condition: 

> \px\, (1.8) 
and hence a natural assumption to be considered concerning the initial conditions p^ and k° is to satisfy 

K°>|p°| on /. (1.9) 

As indicated above, problem (1.3), (1.4) and (1.7), in the case t = 0, has been studied in [15], where a result 
of existence and uniqueness is given using the viscosity /entropy solution framework. Let us just mention 
that in this situation, system (1.3) becomes decoupled and easier to be handled. 

1.3 Statement of the main result 

Remark that the first equation of system (1.3) can be formally rewritten ki = pipx/nx which shows the 
singularity as Kx goes to zero. Nevertheless, for this system we have the following result. 

Theorem 1.1 (Global existence of a solution). Let p^,kP G C°°{I) satisfying (1.9), 

«;°(±1)=±1, p°(±l) = 0, (1.10) 

and the additional conditions: 

DIp\±1)=DIk\±1)=Q, s = 1,2. (1.11) 
Then there exists (p, k) such that for every T > 0; 

(p,k) e (C(7 X [0,T]))2 and p G C\l x {0,T)), 

is a solution of (1-3), (1-4) and (1.7). Moreover, this solution satisfies (1.8) in the distributional sense, i.e. 

Hx>\px\ in V'{Ix{0,T)). (1.12) 
However, the solution has to be interpreted in the following sense: 
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1. K is a viscosity solution of KtKx = ptpx in It '■= I >^ i^^T), 

2. p is a distributional solution of pt = pxx — tKx in It, 

3. the initial and boundary conditions are satisfied pointwisely. 

Remark 1.2 The C°° regularity of and , together with (1-11) seems to be essentially technical and are 
related to an existence result for a regularized system (see Section 3, Theorem 3.1), that we use to prove 
Theorem 1.1. 

1.4 Organization of the paper 

This paper is organized as follows: in Section 2, we present the strategy of the proof. In Section 3, we present 
the tools needed throughout this work. This includes some miscellaneous results for parabolic equations, a 
brief recall to the definition and the stability result of viscosity solutions, and a brief recall to Orlicz spaces. 
In Section 4, we show how to choose the regularized solution. An entropy inequality used to determine 
some uniform bounds on the regularized solution is presented in Section 5. Further uniform bounds and 
convergence arguments are done in Section 6. Section 7 is devoted to the proof of our main result: Theorem 
1.1. In Section 8, some numerical simulations related to our physical model are presented. Finally, Section 
9 is an appendix where we show the proofs of some technical results. 



2 Strategy of the proof 

The main difficulty we have to face is to work with the equation 

i^ti^x = Ptpx- (2.1) 

Since p solves itself a parabolic equation (see the second equation of (1-3)), we expect enough regularity 
on p (indeed p is C^), and then we need a framework where the equation involving k is stable under 
approximation. This property is naturally satisfied in the framework of viscosity solutions (see for instance 
[2] and the references therein). Then, assuming Kx > 0, we interpret k as the viscosity solution of (2.1). 
Assuming (1.9), we wifi indeed show that 

M ■.= Kx- \px\ > 0. 

This is formafiy true because M satisfies: 

Mt = bMx + cM, 

with 

, . s, PxPxx pIx PxxxSgn{px) 

b = Tsgn{px) 2~' ^=~2 ' 

l^x l^x ^a; 

where for suitable boundary conditions, we can (again formally) see that M > 0. In order to justify the 
computations on M, we modify the system and we consider the following paraboHc regularization for e > 
small enough: 

Kl = eKlx + ^^^-rpl in /x(0,oo) 

'^x [2. 2) 

= (1 + £)pL - in /x(0,oo), 

with the initial conditions: 

K,%x, 0) = K°'%x), p%x, 0) = p°''{x), (2.3) 
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where k"''^ and p^'*^ are some regularizations of and p° respectively, and the same boundary conditions: 

.) = x for x = ±1 
p^ix,.) = for a; = ±1. ^ ' ' 

The system (2.2) formally reduces to (1.3) for £ = 0. In fact, system (2.2), (2.3) and (2.4) has (under some 
conditions on the initial and boundary data) a unique smooth global solution (see [16, Theorem 1.1]) for 
ae (0,1): 

(p^K^) e C3+«'^(/ X [0,oo))nC~(/x (0,oo)). 

This result will be recalled in the forthcoming section (see Section 3, Theorem 3.1). The next step is to find 
some uniform bounds (independent of e) on this solution; this is done in particular via: 

(1) an entropy inequality shown to be vahd for our regularized model (2.2); 

(2) a bound on — eK%^ uniformly in e. 

In fact, (1) guarantees the global uniform-in-time control of the modulus of continuity in space of our 
regularized solution, while (2) guarantees the local uniform-in-space control of the modulus of continuity in 
time. The entropy inequality can be easily understood. For instance, for £ = and r = 0, we can formally 
check that the entropy of the dislocation densities 

defined by: 



ji _|_ 



satisfies: 

dS{t) _ net-6-f ^ 
dt J I e+ + 6- - ' 

therefore we get S{t) < 5(0) which controls the entropy uniformly in time. Finally, we need to pass to the 
hmit £ ^ after multiplying the first equation of (2.2) by Having enough control on the regularized 
solutions, we can find a solution of the limit equation using in particular the stability of viscosity solutions 
of Hamilton- Jacobi equations. However, the passage to the limit in the second equation of (2.2) is done in 
the distributional sense. 



3 Tools: miscellaneous results on parabolic equations, viscosity so- 
lution, and Orlicz spaces 

3.1 Miscellaneous results on parabolic equations 

We first fix some notations. Denote 

/t := / X (0, T), 7t := / X [0, T], d^Ir := / U (97 x [0, T]), and || • |Up(o) = || • 
Define the parabolic Sobolev space Wp'^{It) , 1 < p < oo by: 

W^'HIt) ■■= {u e LP(It); {ut,u,,u^^) e (^^'(/t))'} • 

We start with a result of global existence and uniqueness of smooth solutions of the regularized system (2.2), 
with the initial and boundary conditions (2.3) and (2.4). 
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Theorem 3.1 (Global existence for the regularized system, [16, Theorem l-l])- Let < a < 1 

and < e < 1. Let fp'^ , kP'^ satisfying: 



e C°°(/), p°'^(±l)=0, and k°''{±1)=±1, (3.1) 

(3.2) 



r {l+e)pl'^=Tnl'- on dl 
\(l + £)«°^^=rp°'^ on dl, 
and 

nl''>\p°''\ on I. (3.3) 

Then there exists a unique global solution 



(p^/^") G C3+"'^(/x [0,oo))nC°°(/x (0,oo)), (3.4) 
of the system (2.2), (2.3) and (2.4). Moreover, this solution satisfies : 

kI>\pI\ on Ix [0,cxd). (3.5) 

Remark 3.2 Conditions (3.2) are natural here. Indeed, the regularity (3.4) of the solution of equation (2.2) 
with boundary conditions (2.3) and (2.4) imply in particular condition (3.2). 

Remark 3.3 (Uniform L°° bound on p'^ and k'^ ). We remark, from the boundary conditions (2.4) and 
from the inequality (3.5), that: 

lb1iL-(7x[o,oo)) < 1 and ||k''||l=»(/x[o,oo)) < 1- (3-6) 

W'e now present two technical lemmas that will be used in the proof of Theorem 1.1. The proofs of these 
lemmas will be given in the Appendix. 

Lemma 3.4 (Control of the modulus of continuity in time uniformly in e). Let p > 3, and 

g Wp'^{It). Suppose furthermore that the sequences 

(M^)e and (.f )s = (w| - ewDe, 

are locally bounded in It uniformly for e G (0, 1). Then for every V CC It, there exist two constants c > 0, 
£o > depending on V , and < /3 < 1 such that for all < s < eq: 

\-'i-^t + h)-u^{^,t)\ ^^^ V(.,t),(x,i + /.)ey. 

Lemma 3.5 (An interior estimate for the heat equation). Let a G C°°{It) H L^{It) satisfying: 

at = axx on It, 
then for any V CC It, an open set, we have: 

< c||a||i,7.j,, Vl<p<oo, 

where c = c{p, V) > is a positive constant. 
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3.2 Viscosity solution: definition and stability result 

Let n C M" be an open domain, and consider the following Hamilton- Jacobi equation: 

H{x, u{x), Du{x), D^u{x)) =0, Va; e Cl, (3.7) 
where if : x M x M" x MJ^^ ^ M is a continuous mapping. 

Definition 3.6 (Viscosity solution of Hamilton- Jacobi equations). A continuous function u : O i— > 
M. is a viscosity sub-solution of (3.7) if for any <p € C^(0;]R) and any local maximum xq & fl of u — (p, one 
has 

H{xo,u{xo),D(j){xo),D^(j)ixo)) < 0. 
Similarly, u is a viscosity super-solution of (3.7), if at any local minimum point Xq G ^ ofu—<j), one has 

H{xo,uixo),D(l>{xo),Dmxo)) > 0. 

Finally, ifu is both a viscosity sub-solution and a viscosity super-solution, then u is called a viscosity solution. 

To get a "non-empty" and useful definition, it is usually assumed that H is elliptic (see [2]). This notion of 
ellipticity will be indirectly used in Section 7. In fact, this definition is used for interpreting solutions of the 
first equation of (1.3) in the viscosity sense. This will be shown in Section 5. To be more precise, in the case 
where O = It, we say that m is a viscosity solution of the Dirichlet problem (3.7) with u = C, & C{d^lT) if: 

(2) u is a viscosity solution of (3.7) in It, 

(3) u = C on d^lT. 

For a better understanding of the viscosity interpretation of boundary conditions of Hamilton- Jacobi equa- 
tions, we refer the reader to [2, Section 4.2]. We now state a stability result for viscosity solutions of 
Hamilton- Jacobi equations. 

Tlieorem 3.7 (Stability of viscosity solutions, [2, Lemma 2.3]). Suppose that, fore > 0, G C{Cl) 
is a viscosity sub-solution (resp. super-solution) of the equation 

H%x,u^,Du^,D'^u^) =Q in fi, 

where [H^)^ is a sequence of continuous functions. If u locally uniformly in O and if H locally 

uniformly m x M x M" x , then u is a viscosity sub-solution (resp. super-solution) of the equation: 

H{x,u,Du,D^u) = Q in f2. 

3.3 Orlicz spaces: definition and properties 

We recall the definition of an Orlicz space and some of its properties (for details see [1]). A real valued 
function ^ : [0, oo) — > M is called a Young function if 

*(t) = / i;{s)ds, 
Jo 

where ip : [0, oo) — > [0, oo) satisfying: 

• ■^(O) = 0, > on (0, oo), V(i) ^ oo as f — > oo; 

• is non-decreasing and right continuous at any point s > 0. 
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Let \1/ be a Young function. The Orlicz class K^{I) is the set of equivalence classes of real- valued measurable 
functions u on I satisfying 



Ix < +00. 



Definition 3.8 (Orlicz spaces). The Orlicz space L\i,{I) is the linear span of K^{I) supplemented with 
the Luxemburg norm 

Ml^H) = inf ik>0; / VP (\^] < l\ , (3.8) 



I V k 

and with this norm, the Orlicz space is a Banach space. 

The function ^ 

= / (j){s)ds, (j){s) = sup t, 

Jo i/'(*)<s 

is called the complementary Young function of An example of such pair of complementary Young functions 
is the following: 

^-(5) = (l + .s)loo(l + s)-s and $(s)=e"-s-l. (3.9) 

We now state a lemma giving two useful properties of Orlicz spaces that will be used in the proof of Lemma 
5.4. 

Lemma 3.9 (Norm control and Holder inequality, [17]). If u & L^{I) for some Young function ^, 
then we have: 

<1 + J^'b{\u{x)\)dx. (3.10) 

Moreover, ifvG L^{I), $ being the complementary Young function of ^ , then we have the following Holder 
inequality: 



uvdx 



<nu\\L,(I)\\v\\L^(iy (3.11) 



4 The regularized problem 



As we have already mentioned, we will use a parabolic regularization of (1.3), and a result of global existence 
of this regularized system from [16] (see Theorem 3.1). In order to use this result, we need to give a special 
attention to the conditions on the initial data of the approximated system and (see (3.1), (3.2) and 
(3.3)). This section aims to show how to choose the suitable initial data p^'^ and k*^'^ in order to benefit 
Theorem 3.1. Let and k° be the functions given in Theorem 1.1. Set 



„o,...^ _ P°(:^l±£lM and K°-(.) = ^^^" 
(l + £)2 ^ ' l + e 



p^'^ix) = ' \\ ; Z and K"'-(ar) = \ \ \ , (4.1) 
with the function defined by: 



The function </) enjoys some properties that are shown in the following lemma. 
Lemma 4.1 (Properties of (f)) 

The function <j) given by (4-2) satisfies the following properties: 
(PI) ct>,4:\ai = 0; 
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(P2) cl> = 1; 

(P3) \(j)'{x)\ <1/\t\ for x€l. 
Proof. Straightforward computations. 



□ 



Form the above lemma, and from the construction of /9°'^ and k^'^ (see (4.1)) together with the properties 
enjoyed by p° and (see (1.10) and (1.11)), we write down some properties of p°'^ and kP'^ . 

Lemma 4.2 (Properties of p^'^ and kP'^) 

The functions p^'^ and k^'^ given by (4-1) , satisfy the following properties: 



(P4) p°'^(±l) = ±1, and K°'^(±l) = 0; 



(P5) {l + e)Kl'^\g^ = Tpl''\g^ 



and (1 + e)p] 



XX \QI 



= TK. 



0,e| 

X \dl' 



(P6) > \pr\ + 



£(i-MI0'l) 



l + e 

Proof. Straightforward computations 



>\P'. 



□ 



Remark 4.3 (The regularized solution {p^,K'^)). Properties (P4)-(P5)-(P6) of Lemma 4-2 implies con- 
dition (3.1)-(3.2)-(3.3) of Theorem 3.1. In this case, call 



the solution of (2.2), (2.3) and (2.4), given in Theorem 3.1, with the initial conditions 

p{x, 0) = p°'^ and k{x, 0) = 

that are given by (4-1) ■ 



(4.3) 



5 Entropy inequality 

Proposition 5.1 (Entropy inequality). Let {p^,K^) be the regular solution given by (4.3). Define 6»±'^ 
by: 



then the quantity S(t) given by: 



satisfies for every t > 0: 



S{t) = j^^e^'%x,t)\oge^'^{x,t)da 



Sit)<S{Q) + ^. 



(5.2) 



(5.3) 



Proof. From (3.5), we know that k% > \p%\, hence 6^'^ > 0, and the term log(6'^'^) is weh defined. Also 
from the regularity (3.4) of the solution (p^, k'^), we know that 9^''^{.,t) € C(/) for all t>0, hence the term 
S{t) is well defined. We derive system (2.2) with respect to x, and we write it in terms of 9^'^, we get: 



r]9-''+se- 



(5.4) 
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We first remark that: 

Since and pf are zeros on dl x [0,oo), then 

9+'^ -9-'% _^ 



Using (5.5), we compute for t >0: 



T e+'^ + set'' = - 1———^ - T e-^' + eO-^^ = O on 5/ x [0, oo). (5.5) 



- j- + + 

where we have integrated by parts in the second fine. By Young's InequaHty, we have: 
and hence 

Moreover, we have from (2.4), that 
and therefore 

5'(0<y. 

Integrating the previous inequality from to t, we get (5.3). □ 
An immediate corollary of Proposition 5.1 is the following: 

Corollary 5.2 (Special control of k%) 

For all t >0, we have: 

Kl{x, t) \og{K%{x, t))dx < 5(0) + ^ + 2, (5.6) 

where S is given by (5.2). 

The proof of Corollary 5.2 comes from the following inequality. 
Lemma 5.3 For every x,y > 0, we have: 

{x + y) \og{x + y) <x log(x) + y log(y) + x log(2) + y. (5.7) 
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Proof. Direct computations. 



□ 



Proof of Corollary 5.2. From (5.1), it follows that 

ki = e+'^ + e-'^ > 0. 

Then we have for f > 0: 



< 6»+'= log(6»+'=) + e-^^ log(6'-'=) + 6»+'= log 2 + e- 
^ 0+^' \og{0+^') + 0-'^ log(0-'")^ + log 2 + 1 



< S{t) + 2. 

Here we have used Lemma 5.3 with x = 0^'^ and y = 0~''^ for the second line, and we have used for the 
third line, the fact that 



I 0^'' = \j^Kl±pl = i[«^(l, .) - K^-l, .)] = 1. 



Using (5.3), the result follows. 

Lemma 5.4 (Control of the modulus of continuity in space) 

Let u € C^{I), Ux > 0, satisfying 

j Ux\0g{Ux) < Cl, 

for some positive constant Ci, then we have for any x, x + h G I with h > 0; 

2(ci + l + log2) 



□ 



|log/i| 



(5.8) 



■ Ux) Ux 



\u{x + h) — u{x)\ < 

Proof. Let x,x + h G I. 

Step 1. {ux € L^{x,x + h) with ^ given in (3.9)) 

We compute 

j-x+h /■ 
j *(Wx) < y^(l+Ux)l0g(l + l 

< j Ux \og{ux) + log 2 < Cl + log 2, 
where we have used (5.7) in the second line inequality. Hence from (3.10), we get 

\\Ux\\L^(x,x+h) < Cl + 1 + log 2. 

Step 2. (Estimating the modulus of continuity) 

It is easy to check that the function 1 lies in i$(x, x+h) for $ given by (3.9), and that ||l||i,4,(x,a:+/i) ^ ~ 
Therefore, by Holder inequality (3.11), we obtain: 



\u{x + h) — u{x)\ 



x-\-h 



Ux-l 



< 2||Ua;||j;,^(a:^a;_l_;j)||l||2;,^(j:_j:_l_/i) < 



2(ci + l + log2) 
\\ogh\ 
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and the result follows. 



□ 



Remark 5.5 As mentioned to us by Jerome Droniou, it is possible to estimate directly the quantity \u{x + 
h) —u{x)\ < yu^h] splitting the integral J^^'^ Ux on the set where Ux is bigger and lower than A, and then 
optimizing on the parameter A. 



6 An interior estimate 

In this section, we give an interior estimate for the term 

A'=p%- tk'. (6.1) 

that win be used in the passage to the hmit as e goes to zero in the regularized system. We start by deriving 
an equation satisfied by A^: 

Al = {l^e)A%,-^-^Al. (6.2) 

We now show an interior U' estimate concerning the term A^ . This estimate gives a control on the local U' 
norm of A^ by its global V- norm over It, and it will be used in the following section. More precisely, we 
have the following lemma. 

Lemma 6.1 (Interior estimate). Let < e < 1 and 1 < p < oo. Then the quantity A^ given by (6.1) 
satisfies: 

Pllp.y <c(||A^|ji,7, + l), (6.3) 
where V is an open subset of It such that V CC It, and c = c{p, V) > is a constant independent of e. 

Proof. Throughout the proof, the term c = c(p, V) > is a positive constant independent of e, and it may 
vary from line to line. A simple computation gives: 

~'^'i^^x — ~'''~3'(Pxx ~ ''''^x) 
l^x l^x 

_ PxPxx I 2 £ 

= -r{i^l-ei^%x)- (6.4) 

Define as the unique solution of 

on a' It, 

where the existence and uniqueness of this equation is a direct consequence of the theory for parabolic 
equations (see for instance [18, Theorem 9.1]) using in particular the fact that e C^(/t). Moreover, from 
the regularity (3.4) of kF , we can deduce that € C°°{It)- Let 




-t{hI - £K%^), with = - 



We calculate: 



= -^[(1 + ^Xxt + < - £((1 + ^Xxxx + <x)] 
— "'''(1 + ^)i'^xxt ~ ^'^xxxx) ~ '^{'^t ~ ^'^xx) 

= (1 + £)A.,-^A|, 
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where for the first two fine, we have used (6.5), and for the last line, we have used (6.4). In this case, we 
obtain: 

of = Al-Al 

= (l + e)AL-^A^-(l + £)^,, + ^A^ 

= (i + £)(^L-^L) = (i + £)aL> 

where for the first line, we have used the equation (6.2). We apply Lemma 3.5 to the function a^, after doing 
parabolic rescaling of the form a^{x, t) = {x, , we get: 

||a1U<c(l + e)i-^||a-||i,,„ 

and since < e < 1, we finally obtain 

From the definition of a^, and the last inequality, we finally deduce that: 

\m,y<c{\\A^\,j, + \m,j,). (6.6) 

In order to complete the proof, we need to control the term in (6.6). We use the equation (6.5) 

satisfied by to obtain: 

ll^^llp,-fT = tW^I ^ ^'^'xxWp^It 
~ ''"ll'^xa; "I" ^ IIpj^t 

< c{\\Kx\\p,It + II'«1Ip,/t)- (6.7) 
The LP estimates for parabolic equations (see [16, Lemma 2.7]) applied to (6.5) gives: 

then (6.7), together with the fact that < < 1 (see Remark 3.3), implies that: 

hence the result follows. □ 

7 Proof of the main theorem 

At this stage, we are ready to present the proof of our main result (Theorem 1.1). This depends essentially 
on the passage to the hmit in the family of solutions (p*^, n'^) of system (2.2). Since k% ^ 0, we multiply the 
first equation of (2.2) by k% and we rewrite system (2.2) in terms of we obtain: 

f «f < = e<«L + Px^l on It 
\pI = ^pIx + AI on It. 

We will pass to the limit in the framework of viscosity solutions for the first equation of (7.1), and in the 
distributional sense for the second equation. We start with the following proposition. 

Proposition 7.1 (Local uniform convergence) 

The sequences {p^)e, {pt:)s, {i^^)s, {A^)s o^f^d converge (up to extraction of a subsequence) locally 

uniformly in It as e goes to zero. 
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Proof. Let V be an open compactly contained subset of It- The constants that will appear in the proof 
are all independent of e. However, they may depend on other fixed parameters including V. The idea is to 
give an e-uniform control of the modulus of continuity in space and in time of the quantities mentioned in 
Proposition (7.1), which gives the local uniform convergence. The e-uniform control on the space modulus 
of continuity will be derived from the Corollary 5.2 and Lemma 5.4, while the £-uniform control on the time 
modulus of continuity will be derived from Lemma 3.4. The proof is divided into five steps. 



Step 1. (Convergence of and A^) 



From (3.5), wc know that 



< 1. We apply the interior L'p, p > 1, estimates for parabolic equations 



(see for instance [19, Theorem 7.13, page 172]) to the term A- satisfying (6.2), we obtain: 

\\A'\\w?A^v)<(^2\mp,v'' (7-2) 

where V is any open subset of It satisfying V CC V CC It- The constant C2 = C2{p,T,V,V') can be 
chosen independent of s first by applying a parabolic rescaling of (6.2), and then using the fact that the 
factor multiplied by A^^ in (6.2) satisfying 1 < 1 + £ < 2. At this point, we apply Lemma 6.1 for A^ on V, 
we get: 

\\A%y, <c,i\\A^,j, + l), (7.3) 
and hence the above two equations (7.2) and (7.3) give: 

<C4(P1|i,/, + l). (7.4) 

We estimate the right hand side of (7.4) in the following way: 

\\A'h,iT = [ \pl-rK^\ 

J It 

J It 

< (2 + t)T, 

where we have used the fact that \p%\ < k% (see (3.5) of Theorem 3.1) in the second line, and the fact that 
< \k^\ < 1 (see Remark 3.3) in the last fine. Therefore, inequality (7.4) implies: 

< C5, l<p<00. (7.5) 

We use the above inequality for p> 3- In this case, the Sobolev embedding in Holder spaces (see [16, Lemma 
2.8]) gives: 

W^'\V)^C^+"'^{V), a = l-3/p 

and hence (7.5) impHes: 

||^1l^.,.,i^(^^<C6, (7.6) 

which guarantees the equicontinuity and the equiboundedness of {A^)s and {A^)^. By the Arzela-AscoH 
Theorem (see for instance [3]), we finally obtain 

A' — >A and A^ — > A^, (7.7) 

up to a subsequence, uniformly on as £ ^ 0. 

Step 2. (Convergence of k^) 
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We control the modulus of continuity of in space and in time, locally uniformly in e. 
Step 2.1. (Control of the modulus of continuity in time) 
The first equation of (7.1) gives: 



and hence, using the fact that 



Kt — EK^x -|- ^ _ -f^-xi 



< 1, together with (7.6), we get: 



(7.; 



oo,V 



Also, by (3.6), we have: 

||k^||cx),v < 1- 

This uniform bound on together with (7.8) permit to use Lemma 3.4 to conclude that 
\K%x,t) - K^{x,t + h)\<C7hf^, {x,t),{x,t + h) &V, 0</3<l, 



(7.9) 



which controls the modulus of continuity of with respect to t uniformly in e. We now move to control the 
modulus of continuity in space. 

Step 2.2 (An e-uniform bound on 5(0)) 

Recall the definition (5.2) of S{t): 

S{t)= jY^e^'%x,t)\oge^'^{x,t)dx, 
Ji j_ 

with 



3±,£ 



^ Px 



Hence 



5(0) = J^ — 



i-x / X i-x 



0,e _ J),e 



log 



i^x Px 



2 V 2 J Ji 2 °V 2 

Using the elementary identities a;logx < x^ and {x ± y)^ < 2{x^ + y^), we compute: 



5(0) , /( 



+ 



'^x Px 



— llA'a; Il2,7 ^ W^^x \\2,I- 



(7.10) 



From (4.1), we know that: 



and 



_0,e| 



Px + £^<J 



{l+ef 



< 



<\pl\ + h 



l + £ 



Using the above two inequalities into (7.10), we deduce that: 

S{0)<2{\\pl\\lj + K\\lr + 2). 
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step 2.3. (Control of the modulus of continuity in space and conclusion) 

We use the uniform bound obtained for S{0) in Step 2.1, together with the special control (5.6) of k% given 
in Corollary 5.2, we get for all < t < T: 



Kl{x,t) log{Kl{x,t))dx < 2{\\p"JIj + \\4\\li + 2) + ^ + 2, 



therefore 



j^Kl{x,t)\og{nl{x,t))dx <cs, yO<t<T. (7.11) 
Inequality (7.11) permits to use Lemma 5.4, hence we obtain: 

\K^{x + h,t)- K%x,t)\ < jj^^, ix,t),{x + h,t) G It, (7.12) 

Inequalities (7.9) and (7.12) give the equicontinuity of the sequence {K^)e on V, and again by the Arzela- 
AscoH Theorem, we get: 

K, (7.13) 

up to a subsequence, uniformly on ^ as £ ^ 0. 
Step 3. (Convergence of p^) 

As in step 2, we control the modulus of continuity of in space and in time, locally uniformly with respect 
to s. 

Step 3.1. (Control of the modulus of continuity in time) 

The second equation of (7.1) gives: 

hence, from (7.6), we deduce that: 

\\Pt - ^PIx\\oo,V < C6, 

and from (3.6), we have: 

||p^||cx>,y < 1. 

The above two inequalities permit to use Lemma 3.4, we finally get: 

\p%x,t)- p^{x,t + h)\<C7h^, {x,t),{x,t + h) €V, 0<(3<1, (7.14) 
which controls the modulus of continuity of p'^ with respect to t uniformly in s. 
Step 3.2. (Control of the modulus of continuity in space and conclusion) 

The control of the space modulus of continuity is based on the following observation. From (3.5), we know 
that IpII < k| on It- Using this inequality, we get, for every {x,t), {x + h,t) € It- 

j-x+h fX+h 

\p'{x + h,t)-p'{x,t)\< \pl{y,t)\dy< / HUy,t)dy<\K'{x + h,t)-K'{x,t)\- 



Inequality (7.12) gives immediately that: 



\p%x + h,t)-p%x,t)\<y^^, {x,t),{x + h,t)elT- (7.15) 
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From (7.14) and (7.15), we deduce that: 

^ P, (7.16) 

up to a subsequence, uniformly on 1^ as e ^ 0. 

Step 4. (Convergence of p% and conclusion) 

In fact, this follows from Step 1, Step 2, and the fact that 

P%=A' + tk'^P^, (7.17) 
uniformly on 1^ as e — > 0. In this case, we also deduce that 

■A. = Px- TK. 

The proof of Proposition 7.1 is done. □ 
We now move to the proof of the main result. 

Proof of Theorem 1.1. We first remark that is a viscosity solution of the first equation of (7.1): 

k^kI - eK%K%x - p%A% = on It- (7-18) 
Remark 7.2 The equation (7.18) can be viewed as the following Hamilton- Jacobi equation of second order: 

H%X,Dk',D'^k') =0, X = {x,t)€lT (7.19) 

with 

Dk^ = {kI,kI) and D^k^ = 

where is the Hamiltonian function given by: 

-.ItxR^x m2x2 — > R 

"^"^ (7.20) 
iX,p,M) ^ H^{X,p,M)=pip2-epiMi^-pl{X)Al{X), 

P = iPi,P2) and M = (Afij)^ ^=1,2. 

From (7.7) and (7.17), we deduce that {H'^)^ converges locally uniformly in It x x M^^^ to the function 
H given by: 

sym .^21) 
{X,p,M) ^ H{X,p,M)=p^p2-px{X)Ax{X). 

This, togotlior with tlio local uniform convergence of k*^ to k (sec 7.13), and the fact that is a viscosity 
solution of (7.18), permit to use the stability of viscosity solutions (see Theorem 3.7), which proves that k is 
a viscosity solution of 

H{X, Dk, D'^k) = KtKx - PxAx = in It- (7.22) 
We now pass to the limit £ — > in the second equation of (7.1), we obtain 

Pt = Ax€C{lT) in D'(/t). (7.23) 

From (7.22) and (7.23), we get: 

1. «; is a viscosity solution of KfKx = PtPx in It; 

2. p is a distributional solution of pt = Pxx — tKx in It- 
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Inequality (1.12) can be easily obtained by testing against a nonnegative function </> G C^{It)- Finally, 
let us show how to retrieve the initial and boundary conditions. Indeed, the local uniform convergence 
(p^, K^) — > (p, k), together with the uniform control of the modulus of continuity of these solutions: 

• with respect to x near dl x [0,T] by (7.12); 

• with respect to t near 7 x = 0}, away from and 1 by (7.9), 

and the fact that hP'^ p^'^ p'^ uniformly in /, show that (p, k) G (C(Jt))^, so the initial and 

boundary conditions are satisfied pointwisely, and the proof of the main result is done. □ 



8 Application: simulations for the evolution of elastoviscoplatic 
materials 

Motivated by the simulation of the elastoviscoplastic behavior that are formulated by the model of Groma, 
Csikor and Zaiser [12], this section is devoted to write down the equations of the displacement vector u 
inside the crystal when a constant exterior shear stress t (in this section we take r > 0) is applied on the 
boundary walls (see Figure 1). Also, at the end of this section, we present some numerical simulations 
revealing the evolution of a crystal. Here, as we have already mentioned in the introduction, we suppose that 
the distribution of dislocations is invariant by translation in the y-direction. We consider a 2-dimensional 
crystal (Figure 1) with the displacement vector: 

u= (ui,W2) : I — ^K^. 

We rename the coordinates {x,y) as {xi,X2), and we call (61,62) the corresponding orthonormal basis. We 
define the total strain by: 



1 du ' 

^iiW) = 2 + ^^■"j) ■^^^^ '^J"* = g^.' = 

This total strain could be decomposed into two parts as follows: 

£,,(«) =4(«)+ 4, (8.1) 
where e'^{u) is the elastic strain and is the plastic strain which is given by: 

4=p4, (8.2) 

with 

in the special case of a single slip system where dislocations move following the Burgers vector h = ei {Sij 
is the usual Kronecker symbol). Here p is the resolved plastic strain. The stress field a inside the crystal is 
given by: 

aij = 2p4(w) + Xdij Yl ^kkiu) , (8.3) 

\fe=l,2 J 

with A, p > are the constants of Lame coefficients of the crystal that are assumed (for simplification) to be 
isotropic. This stress field a has to satisfy the equation of elasticity: 

j=l,2 ^ 
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Finally, the functions p,K, (solutions of (1.3)) and u are solutions of the following coupled system: 



^ dcFij 



= 



K Pt 



^kk (■") 

\fc=l,2 

\ {djUi + diUj) - pe% 

PtPx 

Pxx '^l^x 



on / X (0,00), 

on / X (0, 00), 

on / X (0,00), 

on / X (0,00), 

on / X (0,00), 



(8.5) 



Equation (8.4) can be reformulated as: 

(A + 2fi)Aui + (A + ^i)di2U2 







(A + 2/i)Au2 + (A + n)d2iui = ndip. 
Here d2P = is due to the homogeneity of the distribution of dislocations in the e2-direction. 



(8.6) 



Calculation of u. We first calculate the value of the displacement u on the boundary walls. Remark first 
that since we are applying a constant shear stress field on the walls, the stress field g there can be evaluated 
as: a ■ n = ±re2, n = ±ei, for x = ±1, 



T 
T 



di. 



Using (8.7) and (8.3), we can derive the following equations on the boundary: 

{diui = on dl, 

fj,{diU2 — p) = T on dl. 

Equation (8.6) leads to the following two equations inside /: 



(8.7) 



(8.8) 



Combining (8.8) and (8.9), with 
finally lead to: 



di[{X + 2ij)diUi]=0 on I 
di{diU2 — p) =0 on I. 

Ml (0,2:2) = U2(0,a;2) = 0, 



(8.9) 



Ui{xi,X2)=0, {xi,X2)€lX 
U2{xi,X2) = —Xi+ f p{x)dx, {xi,X2) & I Xl 

P Jo 



(8.10) 



Formal computation of the long time solution of system (1.3). All computations that will be done 
here are formal. We seek to calculate long time solutions for system (1.3). For this reason, wc first calculate 
the long time solution (stationary solution, nf = pf = 0) of the regularized e-system (2.2), and then we pass 
to the Hmit e — > 0. Doing some computations, the long time {t 00) solution for the £-system is given by: 

p%x,t) = B ^cosh^-:^^^ - cosh ^Y^y^ and K^{x,t) = B sinh ^y^^ (8.11) 
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with B = 1/ sinli y-^j. Passing (again formally) to the limit as £ ^ in (8.11), and using (8.10), we can 
compute the long time displacement U2 inside the material. In fact, we have: 



U2{xi,X2) 



coshr 
sinhr 



Xi + 



sinh TXi 
T sinh T 



(8.12) 



Numerical simulations. The displacement u is numerically computed by discretizing system (1.3) using 
an upwind numerical scheme. The space and time steps in this scheme are well chosen in order to satisfy a 
CFL condition (for the details, see [14, Appendix]). In Figure 2, we show successively the initial state of the 
crystal at time t = Q without any applied stress, then the instantaneous (elastic) deformation of the crystal 
when we apply the shear stress r > at time f = 0+. The deformation of the crystal evolves in time and 
finally converges numerically to some particular deformation which is shown on the last picture after a very 
long time. This kind of behavior is called elasto-visco-plasticity in mechanics. Moreover, on the last picture, 
we observe the presence of boundary layer deformations. This effect is directly related to the introduction 
of the back stress Th = g+j^gt in the model. 




a) i = 0- 



c) t = +00 



Figure 2: Deformation of a slab for model (1.3). 



Remark 8.1 In the last image of Figure 2, it seems that the long time displacement U2 inside the material 
is coherent with our formal computation for the term U2 (see equation (8.12)). 



9 Appendix 

Al. Proof of Lemma 3.4 (control of the modulus of continuity in time) 

Let y be a compactly contained subset of It- Throughout the proof, the constant c may take several values 
but only depending on V. Since V CC It, then there is a rectangular cube of the form 

Q = {XI,X2) X {tl,t2), 

such that V CC Q CC It- In this case, there exists a constant eo, also depending on V such that for any 

< £ < £0, 

and any {x,t) G V, we have: 

(x - 2^/£, X + 2V£) X {t} c Q. 
Moreover, for any (x, t), (x, t + h) ^ V, wc can always find two intervals X and J such that 

{t,t + h) Cl CC J, with {x} X I c Q and {x} x J c Q. 
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Let us indicate that these intervals might have different lengths depending on h and V but we always have 

\J\,\I\ < \t2-tl\. 

Consider the following rescaling of the function defined by: 

u^ix^t) ^u''{^/ex,t). (9.1) 

This function satisfies 

= -"L + {x, t) e (0, X (0, T), 
where f^{x,t) = f^{^/ex,t). Take {xQ,to), {xo,to + h) in V, and let 

Qi = {xa - Ve, .To + Ve) x J and Q2 = (-^0 - 2Ve, + x J. 

These two cyHnders are transformed by the above rescahng into 

We apply the interior , p > 3, estimates for paraboHc equations (see for instance [19, Theorem 7.13, page 
172]) to the function over the domains Qi dd Q2, we get 

\\nwS^\Q.)<<mr>,Q. + yX,Qj- (9-2) 
Using the local e-uniform boundedness of {u^)£, and (/^)£, we get: 

hence, inequality (9.2) implies: 

\\'j'^\\w^•\Q^) ^ c. (9.3) 
We use the Sobolev embedding in Holder spaces (see for instance [16, Lemma 2.8]): 

W^'\Qi)^C'+"'^{Qi), p>3, a=l-3/p, 

to obtain, from (9.3), that: 
and hence 

\u^{xo/y/s,to + h) - U^{XQ/y/s,tQ)\ 

h 2 

then from (9.1), 

|M^(a:o,^o + h) - ^"(3:0,^0)1 

7^+^ — 
h 2 

Choosing /3 = we get the desired result. □ 
A2. Proof of Lemma 3.5 (An interior estimate for the heat equation) 

Recall that a is a solution of the heat equation on It, 
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The proof of lemma 3.5 is a direct computation using a mean value formula for solutions of the heat equations. 
Usually, basic mean value formulae of the solution of the heat equation are expressed through unbounded 
kernels (see for example [7, Theorem 1]), where a can be expressed as: 

a(xo,to) = (47rr2)-i/2 / a{x,t)^^^^^dxdt. (9.4) 

Jnrixo,to) 4(*0 - *) 

Here, {xo,to) e It, {x,t) e ^^(-^^0)1 and r > small enough in order to ensure that the parabolic ball of 
radius r: ^ 

Qrixo, to) = l^{x,ty, to- <t <to, {x-xof <2{to-t)log(^j^^^ cIt- (9.5) 

In our case, wc need a mean value formula similar to (9.4) but with a bounded kernel on flrixa, to). In [10], 
the authors have given such a representation formula for the solution of the heat equation. The following is 
a direct corollary of [10, Theorem 3.1]): 

Corollary 9.1 (Mean value formula with bounded kernels, [10, Theorem 3.1]) 

Let u S C'^(T>) be a solution of the heat equation: 

ut = Uxx on V, 

where V is an open subset o/M^ containing the modified unit parabolic ball fl[.{xo,to), r > 0, with 
fl'.,{xo,to) = ^{x,t); to-r^ <t< to, |a; -a;o|^ < 8(to -*)log 

Then we have: 



to-t 



u{xo, to) = ^ ^1 / u{x, t)E ( — - — , — — ) dxdt, (9.6) 



X — Xo t — to 
|0;(a;o,to)| yn;(xo,to) V r ' r 
where c> 0, \Cl[.{xQ,to)\ = cr^ , and the kernel E satisfies: 



\\E{x,t)\\ < c, (9.7) 

and c > is a fixed positive constant. 

Remark 9.2 In the above corollary, which is an application of [10, Theorem 3.1] in the case m = 3, an 
explicit expression of E can be given by: 



E{x,t) = ^{-x'' + mog{-t)y 



x^ 3(-a;^ + 8tlog(-t)) 
4t2 + 20*2 



where is the volume of the unit ball in . For a more general expression of E, we send the reader to 
]10, Equality (3.6) of Theorem 3.1]. 
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